Anisotropic flow in event-by-event ideal hydrodynamic simulations of ^Snn~- 

Au+Au collisions 
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We simulate top-energy Au+Au collisions using ideal hydrodynamics in order to make the first 
comparison to the complete set of mid-rapidity flow measurements made by the PHENIX Col- 
laboration. A simultaneous calculation of V2, V3, va, and the first event-by-event calculation of 
quadrangular flow defined with respect to the V2 event plane (uffvl^}) gives good agreement with 
measured values, including the dependence on both transverse momentum and centrality. This pro- 
vides confirmation that the collision system is indeed well described as a quark-gluon plasma with 
an extremely small viscosity, and that correlations are dominantly generated from collective effects. 
In addition we present a prediction for vs. 

PACS numbers: 25.75.Ld, 24.10.Nz 



I. INTRODUCTION 

Evidence suggests that in a collision between ultra- 
relativistic heavy nuclei, a strongly-interacting, low- 
viscosity fluid — the quark-gluon plasma (QGP) — is 
created. The clearest indication of this behavior is seen 
in the azimuthal anisotropy 1] among the bulk of emit- 
ted particles. In theory, one characterizes this anisotropy 
in terms of a single-particle probability distribution for 
each collision event. By writing this distribution as a 
Fourier series with respect to the azimuthal angle of out- 
going particles <f>, one can define flow coefficients v n and 
event plane angles ^> n : 
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where the brackets indicate an average over the single 
particle probability, and the event plane angles ty n are 
chosen such that v n are the (positive) magnitudes of the 
complex Fourier coefficients. 

Experimentally, one measures the azimuthal depen- 
dence of event-averaged correlations between detected 
particles. These measurements indicate thepresence of a 
very large "elliptic flow" coefficient w 2 @, 3 , which typ- 
ically can only be reproduced in calculations where the 
system is modeled as a strongly-interacting fluid. In this 
picture, the large momentum anisotropy is generated as 
a hydrodynamic response to the spatial anisotropy of the 
nuclear overlap region in collisions of non-zero impact 
parameter. It even appears that the created quark-gluon 
plasma must be an almost perfect (zero viscosity) fluid, 
with a ratio of shear viscosity to entropy density r\j s that 
is at most a few times l/(47r), a value that was famously 
conjectured to be a universal lower bound [4]^j. However, 



the extraction of a precise finite value is hampered by 
poor knowledge of the earliest stages of the collision, as 
well as other uncertainties [10| . 

An important recent development was the realization 
of the importance of quantum fluctuations, which in par- 
ticular implies an event- by- event breaking of the sym- 
metry naively implied by the collision of identical nuclei. 
Specifically, the coefficients v n are generally non-zero also 
for odd n [ll|, the event plane angles do not necessar- 
ily poin t in the same direction as the impact parame- 
ter [T2I, LL3|; an d these quantities fluctuate significantly 
from one event to another, even at a fixed impact pa- 
rameter pL4j ] . 

These insights led to the possibility that all of the 
measured long-range correlations may be generated solely 
from collective behavior pj], [l5| . 

Several new flow observables — specifically ones im- 
plied by the presence of event-by-event fluctuations - 
were recently measured for the first time |16l - [2fj |. Stud- 
ies of these new observables indicate that, individually, 
they indeed appear to have properties that are consis- 
tent with a hydrodynamic origin [15|, I21H23J . However, 
they have not yet all been reproduced in a single calcu- 
lation within one model using a single set of parameters. 
This has left some lingering doubt about whether the 
interpretation in terms of collective behavior is indeed 
correct [24|. In addition, each measurement provides an 
independent constraint on theory, so identifying models 
and sets of parameters that can simultaneously satisfy 
all the constraints is a necessary first step in reducing 
various theoretical uncertainties. 

In this Letter we perform state-of-the-art ideal hydro- 
dynamic calculations and compare the results to the first 
measurements [16[ of these new observables at Relativis- 
tic Heavy-Ion Collider (RHIC) as well as previous mea- 
surements by the same collaboration 25]. Other groups 



1 The bound is now known to be violated in some theories [5|4j|, 



and it may even be possible to have an arbitrarily small value 
though the effective viscosity may still have a finite bound 



have presented calculations from some of these observ- 
ables using event-by-event ideal [26J - 4291] or viscous |30| 
hydrodynamics, or transport models [31j , The present 
study encompasses simultaneously, for the first time, all 
the measured flow observables at midrapidity. 



II. OBSERVABLES 

All the experimental results considered here were ob- 
tained using the event-plane method [32j. With this 
method, one first identifies an event plane \l/„ in each 
event using a specific detector at forward rapidity, and 
then calculates the correlation of particles near midra- 
pidity with this event plane, e.g., 



v n {* n } = (cosn(<f> -*„)), 



(3) 



where the brackets indicate an average over particles in 
a large number of events. A rapidity gap with the event- 
plane detector suppresses nonflow correlations [15|, I33J . 
At RHIC, "triangular flow" wsj^a} and "quadrangu- 
lar flow" ^{vE^} were measured for the first time, as a 
function of the particle transverse momentum p t in vari- 
ous centrality classes by the PHENIX collaboration [16| 
(preliminary data from STAR have now also been pre- 
sented [H). 

Previously, a different quadrangular flow observable 
has been measured, defined with respect to ^2 |25l.[34|. 
We use a different notation for this quantity to avoid 
confusion: 



u 4 {* 2 }= (cos4(0-# 2 ))- 



(4) 



v n is analyzed using a large sample of events, and its 
value fluctuates from one event to the other. These fluc- 
tuations (which were not appreciated when the method 
was developed) , combined with the use of a finite sample 
of particles in the analysis, cause the measured value to 
deviate from the event average of the theoretical coeffi- 
cients defined in Eq. (p}. Generally, v n {$> n } lies between 
the mean value and the root-mean-square (rms) value 
of Vn- One can parameterize the resulting measurement 
as 1351: 
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where here the brackets indicate an average over events. 
The value of a de pen ds on the event plane resolution 
Resj^n} ~ v n y/N [36}: If the resolution is poor, a ~ 2, 
and the measured v n is a rms value, while if the resolution 
is large, a ~ I, and the result gets closer to the mean 
value. 

The most recent data from PHENIX has a maximum 
event plane resolution of 0.74 (for t> 2 around 30% cen- 
trality 25]) and much smaller for w 3 and V4 [16j . which 
implies a > 1.81 [3a|- So in general the results are very 
close to a rms value of v n . Nevertheless, in the follow- 
ing we compute both limiting cases a = 2 and a — 1 



in order to show the size of the effect of fluctuations on 
event-plane analyses. 

Likewise, the measured value ^4(^2} depends on the 
resolution [37| . and is usually close to (v±v\ cos(4* 4 



^2))/vK)i but with increasing resolution approaches 

(w4COs(4*4-4^2))- 



III. RESULTS 

Using the hydrodynamic code NeXSPheRIO [38|], we 
simulate top-energy Au-Au collisions at RHIC. This code 
solves the equations of ideal relativistic hydrodynamics 
using fluctuating initial conditions from the event gener- 
ator NeXus J39J]. 

NeXus aims at a realistic and consistent approach of 
the initial stage of nuclear collisions [39] . It is a Monte- 
Carlo generator which takes into account not only the 
fluctuations of nucleon positions within nuclei J30j . but 
also fluctuations at the partonic level: the momentum of 
each nucleon is shared between one or several "partici- 
pants" and a "remnant" , which implies non-trivial dy- 
namical fluctuations in each nucleon-nucleon collisions. 
The resulting full energy-momentum tensor is matched 
to a hydrodynamic form, resulting in a fluctuating flow 
field in addition to a fluctuating initial energy density, in 
all three spatial dimensions, with the transverse length 
scale of the fluctuations set mostly by the size of the in- 
coming nucleons. 

At the end of the hydrodynamic evolution, discrete 
particles are emitted using a Monte-Carlo generatoiQ 
NeXSPheRIO provides a good description of rapidity and 
transverse momentum spectra [431 ] . elliptic flow v 2 [44|, 
and the rapidity-even v\ observable, directed flow at 
midrapidity |45| . In addition, it is known to reproduce 
the lon g-ra nge structures observed in two-particle corre- 
lations [46]. All parameters were fixed from these earlier 
investigations, before any of the new observables (^3, V4) 
were measured — nothing has been tuned here. 

For this work, we generated 110 NeXus events each in 
5% centrality classes up to 60% centrality, solving the hy- 
drodynamic equations independently for each event. As 
in Ref. 47J, at the end of each hydro event, we run the 
Monte-Carlo generator many times, so that we can do 
the flow analysis using approximately 6 x 10 5 particles 
per event. This significantly reduces statistical noise and 
allows for an accurate determination of v n and \& n in ev- 
ery event. It also suppresses non-flow correlations from, 
e.g., particle decays. These quantities are then calculated 
by Eq. ([3]), with the average taken over all particles in 
the pseudorapidity interval — 1 < r\ < 1. The procedure 
used to measure v n in hydrodynamics thus mimics the 



2 Freeze-out occurs at a constant temperature. Hadrons do not 
rescatter after freeze-out I40I442H . but resonance decays are im- 
plemented. 
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FIG. 1. (Color online) Results for V n {^n} for n — 2-5, compared to published data from the PHENIX collaboration [la ]. 
Closed and open symbols correspond to two different ways of averaging over events (mean and rms value, respectively). Error 
bars represent statistical uncertainty from the finite number of events. The left column (0-10%) represents the 10% most 
central collisions, which each column to the right increasingly peripheral. 



experimental procedure, with two differences: (i) there is 
no need for a rapidity gap, because nonflow correlations 
are negligible; (ii) there is no need for a resolution cor- 
rection, because the huge multiplicity per event ensures 
that the resolution is close to 1 for all events |47| . 

Fig. Q] displays v n as a function of the particle trans- 
verse momentum p t , averaged over events in a centrality 
class. The average over events is estimated in two differ- 
ent ways in order to illustrate the effect of event-by-event 
flow fluctuations on the experimental analysis. The first 
estimate, labeled NeXSPheRIO-, is a plain mean value 
(corresponding to a = 1 in Eq. ©). The second esti- 
mate, labeled NeXSPheRIO+ is a weigthed average 
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(6) 



The average of v^l'Fn} over p t is the rms average o f v n 
(a = 2 in Eq. ([5])). For Gaussian flow fluctuations [48], 
the ratio of the rms to the mean is y&pK ~ 1.13 for v% 
and W5, and closer to 1 for V2 and U4. 

Fig. Q] shows that our event-by-event ideal hydrody- 
namic calculation reproduces well the observed central- 
ity and transverse momentum dependence of vi , vj, and 



V4. The pt dependence is a generic feature of ideal hy- 
drodynamics [2l|. The magnitude and centrality depen- 
dence of v n , on the other hand, depend on the initial 
conditions: vi is mostly driven by the almond shape of 
the overlap area, which depends on the particular model 
used [49(|, while hi gher harmonics are mostly driven by 
initial fluctuations [ll| . which explains why they have a 
mild centrality dependence [501 ] - 



Originally, quadrangular flow W4 had been measured 
with respect to the event-plane of elliptic flow. Recent 
results show that ^{"J^} [25] is smaller than 1/4(^4} 16], 
typically by a factor 2 for peripheral collisions, and by 
a factor 5 for central collisions. Although W4JV1/2} is 
smaller, it is measured with a better relative accuracy 
than V4{^ 4 }, because of the better resolution on \E , 2- 
This makes vA^<i\ a useful quantity for detailed model 
comparisons J51J. As in the case of Vnl^n}, we per- 
form the average over events in two different ways in 
order to illustrate the effect of event-by-event flow fluc- 
tuations. The first estimate, labeled NeXSPheRIO-, is 
a plain mean value, as in Eq. (U). The second estimate, 
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FIG. 2. (Color online) Results for the first cvent-by-event hydrodynamic calculation of «4{5 , 2}/i , 2{^ , 2} 2 , compared to published 
data from the PHENIX collaboration [2g]. As in Fig.[T] closed and open symbols correspond to two different ways of averaging 
over events (see text), error bars represent statistical uncertainty from the finite number of events, and smaller percentile refers 
to more central collisions 



labeled NeXSPheRIO+ is a weigthed average: 
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The actual event-plane value is expected to lie between 
these two limits, depending on the resolution [371 ]. 

Since ^{^j can be generated by elliptic flow as a 
second order effect [52|, we scale it by i^l 1 !^} 2 for each 
p t . Hereafter, we denote ujj'I^} and «2{*2} simply by 
v 4 and v 2 - Fig. [2] displays this first event- by- event hy- 
drodynamic calculation of v^/v\ as a function of pt for 
different ccntralities. The measured ratio is remarkably 
constant as a function of pt , and increases mildly for cen- 
tral collisions. Ideal hydrodynamics predicts v^/v\ ~ 1/2 
at high p t for a single event [52]. For all centralities, the 
measured value of v^/v\ is greater than 1/2, even at high 
pt- This can be explained [37j by V2 fluctuations, ex- 
cept for the two most central bins, where one expects 
Vi/v% ~ 1 [371 ] . smaller than the measured value, which 
is between 1.5 and 2 for the most central bin. For these 
two central bins, our results from event-by-event hydro- 
dynamics are in good agreement with experiment (first 
two panels in Fig. (2). This shows that other sources of 
flow fluctuations, other than v-i fluctuations alone, con- 
tribute to i>4/ff . A similar finding has been reported in 
a transport calculation with v^ and i>2 both defined with 
respect to the direction of the impact parameter [3lj j . 
Our calculated W4/W2 is slightly higher than data for the 
next two bins (10 — 20%). Above 20% centrality, data 
are within the range spanned by our calculations. 

The calculations shown here simulate the system evolu- 



tion using ideal hydrodynamics, i.e., with negligible vis- 
cosity. These results prove that no non-zero QGP vis- 
cosity is required to reproduce these data. In fact, this 
calculation requires a negligible viscosity — keeping ev- 
erything else fixed, a viscosity the size of the conjectured 
bound rj/s = 1/Att would suppress v n and destroy the 
remarkable fit to data. In addition, the ratio va/v^pt) 
depends strongly on rj/s, and any non-zero value usually 
tends to destroy the flat curve that ideal hydrodynamics 
predicts [5ll [531] . 

However, this requirement of negligible viscosity de- 
pends crucially on aspects of the model which are not 
entirely constrained. In particular, although the NeXus 
model provides an honest effort at a reasonable descrip- 
tion of the physics, with many realistic elements, there 
is considerable uncertainty about the early stages of a 
heavy-ion collision and the resulting initial conditions for 
hydrodynamic evolution. In principle, another model, 
coupled to viscous hydrodynamics might well be able to 
fit these data. For example Ref. J30J] presents event-by- 
event viscous hydrodynamic calculations with Glauber 
initial conditions that require a value close to rj/s — 0.08 
to give reasonable agreement with the quantities in Fig.Q] 
at several centralities, though they underpredict 1)3 for 
central collisions. Secondly, although v^/v\ is not very 
sensitive to the initial conditions, the effect of non-zero 
viscosity depends significantly on the way it is imple- 
mented at freeze out [531 ] , and the correct implementation 
is an open issue. 

Thus, this work is only a first step in identifying mod- 
els that are compatible with data, and strong conclusions 



cannot yet be drawn about, e.g., the precise value of rj/s. 
Although the success of these calculations are an impor- 
tant milestone, proving that at this point no lower bound 
can yet be placed on rj/s, we can not yet make a precise 
statement about an upper bound - only that it still ap- 
pears unlikely that a value significantly larger than 1/47T 
will be possible. 



IV. CONCLUSIONS 



in a single calculation, this provides a complete, unified 
picture of the bulk evolution of a heavy-ion collision as 
an extremely low-viscosity fluid. Indeed, for our model 
of initial conditions, a negligible viscosity is required for 
a good fit to all mid-rapidity flow observables. Therefore 
no lower bound can currently be placed on the shear vis- 
cosity of the quark-gluon plasma. Further study will be 
needed to determine a reliable upper bound, but finding 
models (such as this one) that are compatible with all 
measured data is a significant first step. 



Using an ideal hydrodynamic model with fluctuating 
initial conditions, we have performed the first simultane- 
ous calculation of V2{ ^2}, ^3(^3}, v^^f^} and uil^} as 
a function of transverse momentum and centrality. Our 
results are in good agreement with the most recent ex- 
perimental results for all the observables at RHIC, at all 
centralities and in a wide range of transverse momen- 
tum. This provides convincing confirmation of the cur- 
rent paradigm that collective effects alone can explain all 
long-range correlations in the soft sector. Further, since 
all such measured correlations are generated consistently 
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